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PACS 05 . 60 . Gg - Quantum transport 
PACS 44 . 10 . +i - Heat conduction 

PACS 66 . 10 . cd - Thermal diffusion and diffusive energy transport 

PACS 85. 75. -d - Magnetoelectronics; spintronics: devices exploiting spin polarized transport or 
integrated magnetic fields 

Abstract. - The heat conduction in an open transverse-field Ising chain is studied by using 
quantization in the Fock space of operators in the weak coupling regimes, i.e. the coupling is 
much smaller than the transverse field. The non-equilibrium steady state is obtained for large 
size systems coupled to Markovian baths at its ends. The ballistic transport is observed in the 
uniform chain and normal diffusion in the random- exchange chain. In addition, the ballistic- 
diffusive transition is found at the intermediate disorder regime. The thermal conductivity k is 



also calculated in the low and high temperature regimes. It is shown that k decays as k • 
high temperatures. 
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Introduction. — Heat and spin transport behaviors 
of one-dimensional systems have intensively been stud- 
ied in both classical and quantum mechanical context for 
several decades [TMl3|. Some classical nonlinear systems 
of interacting particles usually show a diffusive behavior, 
which satisfy Fourier's law, J = — kVT (k is the thermal 
conductivity), relating the macroscopic heat flux to the 
temperature gradient |14H16| . In the quantum regime, it 
is also an important issue that how normal energy or heat 
, is transported through a sample on a microscopic level. 
The ballistic behavior has been proposed in the integrable 
quantum systems, implying that the current-current cor- 
relation functions typically do not decay to zero [7]. The 
chaotic dynamics of a nonintegrable system may yield a 
normal diffusion. One natural question is how the heat 
transports through a quantum disordered system. 

There exists two general theoretical approaches for a 
description of non-equilibrium open quantum spin sys- 
tems. One is the non-equilibrium green's function method 
[T71IT8] . The other is the quantum master equation [I9|l20] . 
In the latter case, various approximations schemes are em- 
ployed, such as Markov approximations, Born approxima- 
tions, secular approximations, and weak coupling approx- 
imations for the system-environment coupling. Some ef- 
fective dissipative equations of motion for reduced density 



matrix p s of the open systems are then derived. 

By using the proper quantum master equation, some 
interesting and important thermal properties have been 
found in the open spin-chain systems by using the Monte 
Carlo wave function method [5"I l2"TH2"3"] . the matrix prod- 
uct operator method [23], and fourth order Runge-Kutta 
method [TJ[31II]. However, the largest system sizes re- 
ported in the literature is smaller than N = 20, and the 
convincing results for the properties and the type of the 
transport arc still lacking to date. In the meantime, some 
progresses in the methodology have also been achieved, 
such as the adaptive time-dependent density matrix renor- 
malization group, the numerically exact diagonalization 
and the quantization in the Fock space of operators etc. 
25 29., which could be applied in the large system suc- 
cessfully. Therefore, the studies on the large size system 
have become an essential and intensive issues. Recently, 
Prosenei al. [33] proposed a method to solve explicitly the 
Lindblad master equation for an arbitrary quadratic sys- 
tem of n fermions in terms of diagonalization of a An x An 
matrix. This method has been successfully applied to the 
far from equilibrium quantum phase transition [30 , 31 in 
one-dimensional XY spin chain for size larger than 100. 

In this paper, we study the heat current properties of the 
well known ID transverse field Ising model with very large 
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sizes. The exchange couplings are considered to be both 
uniform and random. The decisive conditions of the bal- 
listic transport and the normal transport are given. The 
paper is organized as follows: In Sec. II, we describe the 
model and the scheme to solve the quantum master equa- 
tion in detail. Then the numerical results are obtained in 
Sec. III. The conclusion is given in the last section. 

MODEL AND METHOD. - The Hamiltonian for 
an open ID transverse field Ising chain reads 



and T2 ,R = 7rAsTii i ij(a;)7(a;)| w= 2/ l = Xn L:R (2h). Here, 
n-L,/?.M = ( e u/TL - R - l)" 1 is the Bose-Einstein distri- 
bution function (kg — 1)> is the system-bath cou- 
pling strength, and I{oS) denotes the spectral density of 
an Ohmic bath that we choose. The Hamiltonian (1) is 
conveniently expressed as a quadratic form H — w- Hw 
in terms of 27V Hermitian Majorana operators 
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where TV is the number of spins, the operators and <r„ 
are the Pauli matrices for the nth spin, J n is the coupling 
parameter between the nearest-neighbor spins, and h is 
the transverse magnetic field. Here we take h = 1. For 
the disordered system, J n is chosen to distribute on a in- 
terval (0.05,0.15) uniformly, modeling the weak coupling 
condition. For the uniform system, we take J n = 0.15. 
Considering the two thermal baths at two ends and the 
coupling with the Ising chain, the total Hamiltonian can 
be written as 



satisfying the anticommutation relation {w n ,w m } = 
25 n _ m . H is an 2N x 27V antisymmetry Hermite matrix 
(H T = H). Based on the previous transformation, we 
can rewrite Hamiltonian (1) in terms of Majorana fermions 
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where ^l(_r) = Ylk ^kO-k^k is the left(right) phonon 
bath with a k (a,k) the phonon creation(annihilation) oper- 
ator, and H int = crf (ir) Hk \9k a k + 9k a \) is tne interac- 
tion between the chain and the baths. If the coupling 
is weak, a quantum master equation for the systems evo- 
lution can be obtained from our microscopic Hamiltonian 

model by using the usual Born-Markov approximations a 2j _i = — (cj+c 1 ), a,2j 
and the secular approximation [19] . x 

The quantum master equation in the weak internal cou- 
pling limit ({ Jn} *C h) reads (we set Ti = 1) 



where W = w\W2 • • • W2N is a Casimir operator which 
commutes with all the elements of the Clifford algebra 
generated by Wj [35]. Note that W 2 = 1, so it does not 
affect the quadratical system. For convenience, we take 
W = 1. 

Then we construct A N dimensional Pauli algebra with 
a Fock space of operators describing 27V adjoint fermions 
(a- fermions), with an orthonormal canonical basis \P a ) = 



), a 3 e {0,1}. With the definition 
Pa) — $a - AwjPa), the quantum Liouvillean (4) be- 
comes bilinear j£> = a ■ Aa + constl in Hermitian maps 
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The 47V x 47V matrix A can be expressed in a block 



dp 
~dt 



-i[H lP } + D L (p) + D R (p), 



(3) 



tridiagonal form in terms of 4 x 4 matrices as 
/ B L - 2hR Ri 
-2KB. 



where the dissipator Dl refers to the left heat bath and 
D R to the right one, depending on the full density operator 
p of the Hamiltonian (1). Eq. (3) can be rewritten as the 
Lindblad master equation 
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where L u s are the Lindblad operators, representing cou- 
plings to different baths. The weak bath coupling is taken 
into account here. The simplest nontrivial bath opera- 
tors acting only on the first and the last spin are chosen 
(p =1,2,3, and 4) 
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have rf ' R 



«„ ± Refer to Refs. (H[29l[30], we 
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Fig. 1: (color online) Energy profiles E n (odd sites) for a ran- 
dom chain with different sizes N. J n £ (0.05,0.15), Tl — 5, 
Tr = 2. For comparison, the energy profile E n of an ordered 
chain with J n = J = O.f is also listed. The inset denotes the 
local energy E2 versus average temperature T for the disor- 
dered system (N = 4) with different temperature difference A. 
All results are averaged over 200 realizations of disorder. 



The eigenvalues of AN x AN antisymmetric A called 
rapidities can be list in the form of the pairs 
-/3i,/3 2 , -&>, • ■ • ,p2N, -P2N, Re/3, > 0. The corre- 
sponding eigenvectors v p ( p = 1, . . . , AN ) can be defined 
by Av 2 j_ 1 = PjV 2 j_i and Au 2j = —/3jV 2 j- v can be nor- 
malized by using u-y-i ■ 2L 2 j = 1 an d U p -v q — otherwise, 
which can be used to calculate any quadratic physical ob- 
servable in the non-equilibrium steady state (NESS). The 
expectation value is given by 



(WjW k ) 
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NUMERICAL RESULTS. - First, to show the 
correctness of this method, we numerically calculate the 
local energy E n . The local energy density operator reads 

H n = iJ n W2nW 2 n+l + i^(w 2n -iW 2n + W 2n+ iW 2n+2 ) , (13) 

and local energy density is defined as E n = (H n )NESS- 
The temperatures of the left and right baths are Tl = 
T(l+A) and T R = T(l- A), where A is the dimensionless 
temperature difference and T is the average temperature. 
We set the coupling constants uniformly distribute on the 
interval J n e (0.05,0.15), which are the same as those in 
Ref. l.i. A is taken as 0.005, which is equivalent to the 
parameter a = 0.01 in the Ref. pp. The numerical results 
for local energy E n with different sizes are presented in 
Fig. 1. Here we only plot the local energy on the odd sites 
for the similar results on the even sites. Note that there is 




Fig. 2: (color online) Heat current vs temperature difference 
for the uniform system. The average temperature T changes 
from f to 3. The other parameters are h = 1, J„ = 0.f5, 
A = 0.005, and N = 100. 



a clear intersection of the energy profiles E n (disordered) 
for different sizes at the central part of the chain. And 
the inset shows E 2 versus T for N — A. The local energy 
decreases with the increase of A. When A = 0, the results 
can be checked with the canonical one pQ . It is interesting 
that the present local energy is nearly the same as those in 
Ref. 1 . It should be pointed out here that the method is 
not suited to the strong coupling case, since the Lindblad 
equation is only valid for small {Jn}- 

Then we evaluate the heat current in the spin chain in 
NESS 

Qn = i[H n ,H n +i] 

= i[-J n hw 2n w 2n+2 - hJ n+ iw 2n+ iw 2n+3 ]. (1A) 

The heat current as a function of the temperature differ- 
ence for N = 100 chain is presented in Fig. 2 for the 
uniform case and Fig. 3 for the disordered case with the 
mean temperature ranging from T = 1 to 3. According 
to Eq. (5), the high bath temperature can enhance the 
coupling strength between the baths and the spin chain, 
which facilitates the heat transfer. So k increases with the 
augmentation of the average temperature T in the given 
parameters regime, as shown in Fig. 2. For the disordered 
case, k changes little in small A regimes (—0.1 < A < 0.1), 
as shown in Fig. 3. From the order of magnitude of the 
heat current, one can see that the current is easily blocked 
for the disordered case. 

We turn to discuss the classification of heat transport 
properties. Note that a finite current within an infinite 
system demonstrates ballistic transport behavior. In the 
previous studies, the ballistic behavior is observed in the 
integrable system, and the diffusive transport occurs for 
the disordered system. But these conclusions were built on 
the numerical simulations on small systems [TJ[3l[5l[T0l[TT] . 

To classify transport properties more convincingly, we 
simulate the system with the size up to N > 100 for both 
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Fig. 3: (color online) Heat current vs temperature difference 
for the disordered system. The range of T is set from 1 to 3. 
The other parameters are h = 1, J n € (0.05, 0.15), A = 0.005, 
and N = 100. The results are averaged over 100 realizations 
of disorder. 



cases. The numerical results are plotted in Figs. 4 and 
5. Note that the ballistic transport occurs at the uniform 
case, as shown in Fig. 4. The heat current increases with 
the bath coupling strength A. Different from the observa- 
tion in a Heisenberg chain [5], the present heat current is 
only sensitive to the parameter A and independent of the 
system size, which is just the characteristics of the ballistic 
transport only emerging in a quantum system. The known 
Fourier's law is obviously invalid. For the disordered case, 
the current is gradually reduced with increasing the sys- 
tem size, demonstrating the diffusive transport, as exhib- 
ited in Fig. 5. The similar results are also obtained in the 
mass-disordered harmonic crystals 32 . For the present 
system size N — 100, l/N scaling behavior of the heat 
current is observed for A = 0.0025, demonstrating that 
the Fourier's law holds in this case. So we provide an 
evidence of the macroscopic heat transport in the quan- 
tum disordered system. With the increase of the coupling 
parameter A, the linear fit gradually deviates from the ori- 
gin of the coordinate, implying that the finite-size effect is 
more obvious for large coupling parameters, thus the data 
for large systems would be essential to get the correct scal- 
ing for strong coupling. 

At this stage, it should be expected that the ballistic- 
diffusive transition might occur at the intermediate dis- 
order regime. The disorder strength is introduced in the 
coupling parameter J n = 0.1 + Qdis£, where e distributes 
uniformly on the interval (—0.5,0.5). We plot the heat 
current vs reciprocal chain length for different disorder 
strengthes in Fig. 6. The linear fit of the data shows that 
the transition occurs at Qdis = 0.04 ~ 0.05. 

Finally, we obtain the fully thermal conductivity k as 
a function of temperature in the different large sizes for 
the disordered case. In the high temperature regime, it is 
observed that K decays slightly faster than T -2 , as shown 
in Fig. 7. k reaches a maximum value around T = 1.0 



Fig. 4: Heat current versus reciprocal chain length for differ- 
ent bath coupling strength A for the uniform system. The 
parameters are J n = 0.15, h = 1, 0l = 1/Tl = 0.25 and 
Pr = 1/Tij = 0.75. 



with different sizes N . The thermal conductivity decreases 
with increasing system size. The inset presents the rela- 
tion ln{T 2 n) oc 1/T in the low temperature regime. The 
thermal conductivity k and the specific heat c v of a single 
spin have a similar temperature dependence. These ob- 
servations are consistent with the analytic and numerical 
results of Ref. 1 based on a smaller system. 

Conclusions. — In this paper, we have studied the 
heat transport behaviors of an open Ising chain within 
master equation formalism by using quantization in the 
Fock space of operators. The classification of the trans- 
port properties is performed in large size system (over 100 
sites). Compared with the Monte Carlo wave- function 
method, the precision has been considerably improved in 
the present approach. We confirm the ballistic current in 
the uniform system which is integrable. The normal trans- 
port is clearly observed in the disordered system. The bulk 
conductivity decreases with the increase of the system size. 
The heat current exhibits a diffusive behavior above a crit- 
ical interaction strength, which follows Fourier's law in the 
normal transport. Moreover, the ballistic-diffusive transi- 
tion occurs at the intermediate disorder regime. It is also 
observed that the thermal conductivity k has the simi- 
lar temperature dependence as the specific heat c v in the 
weak coupling regime. 
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Fig. 5: (color online) Heat current versus reciprocal chain 
length for the disordered system. System parameters: J n € 
(0.05, O.f 5), h = 1, p L = 1/Tl = 0.25 and p R = 1/T R = 0.75 
. The sizes of the system range from 60 to 120 every other 
10 sites. The results are averaged over 2000 realizations of 
disorder. 
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Fig. 7: Thermal conductivity k versus T for the disordered 
system with J n £ (0.05,0.15), h = 1 and different sizes N = 
60, 80, 100. The bath temperatures are T L = T(l + A) and 
Tr = T(l - A) (A = 0.02). The results are averaged over 200 
realizations of disorder. The inset shows the low temperature 
property of thermal conductivity. 
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